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Robust porous silicon substrates were employed for generating interconnected networks of super- 
conducting ultrathin Nb nanowires. Scanning electron microscopy analysis was performed to inves- 
tigate the morphology of the samples, which constitute of polycrystalline single wires with grain size 
of about 10 nm. The samples exhibit nonzero resistance over a broad temperature range below the 
critical temperature, fingerprint of phase slippage processes. The transport data are satisfactory 
reproduced by models describing both thermal and quantum fluctuations of the superconducting 
order parameter in thin homogeneous superconducting wires. 



PACS numbers: 74.78.Na, 73.63.Nm 



The research activity in the field of superconductiv- 
ity at reduced dimension has been continuously growing 
[JlH, due to the major implications that these results can 
have for understanding fundamental phenomena 0, 3| as 
well as forpossible applications in superconducting elec- 
tronics 043 ■ The behavior of the so-called superconduct- 
ing nanowires, namely ultrathin samples with dimensions 
comparable with the superconducting coherence length, 
£, is dominated by both thermal activated phase slip s 
(TAPS) and quantum phase slips (QPS) processes flOj . 
causing the wires to remain resistive much below the su- 
perconducting transition temperature. The most chal- 
lenging aspect of the experimental study of nanowires is 
the difficulty of fabricating homogeneous samples, since 
it has been widely demonstrated that sample inhomo- 
geneity can be the source of broadened superconduct- 
ing transitions [IlT - fl3l |. From the early studies realized 
on crystalline superconducting whiskers [lij , advances in 
nanofabrication techniques have allowed the realization 
of high quality single crystals nanowires [l5j . Even more 
sophisticated is the approach of using suspended carbon 
nanotubes or DNA molecules as templates for the forma- 
tion of superconducting nanowires |16l4l8l |. A radically 
different approach to nanostructures fabrication based 
on self-assembled growth attracted much attention also 
in the superconducting nanowire field [l9 22]. Most of 
these works rely on the deposition within the channels 
and cavities of porous membranes. Self-assembled meth- 
ods are versatile and reliable bottom-up techniques for 
generating low-cost patterns of nanostructures, assuring 
a highly reproducible geometry on very large areas. 

In this work the formation of interconnected networks 
consisting of Nb ultrathin superconducting nanowires is 
achieved by using Porous Silicon (PS) as template sub- 
strate. The extremely reduced film thickness favors the 
deposited material to occupy only the substrate pitch, 
therefore the as-sputtered films result as a network of 
interconnected wires, whose average width, w, can be as- 
sumed to be equal to the periodic pore spacing minus the 
pore diameter. Due to the extremely reduced character- 



istic dimension of the substrate w is comparable to the 
superconducting coherence length £, and hence each in- 
dividual wire behaves as a one-dimensional (ID) object. 
As a consequence, the whole samples show broadened 
resistive transitions, which can be described by theoreti- 
cal models for both thermal and quantum fluctuations of 
the order parameter in ID superconductors 23 25] . The 
proposed technique allows one to overcome the problem 
of handling f rag ile membranes usually used at this pur- 
pose [ill [20, [Hj employing robust PS crystals as stable 
support. Moreover, the analyzed systems are rather sim- 
ple and macroscopically large objects, which, however, 
reveal fascinating quantum effects. Finally, compared to 
similar self-assembled nanowire networks [22J exhibiting 
thermal phase slippage, the system presented in this work 
provides strong evidence of QPS. 

PS was fabricated by electrochemical anodic etching of 
n-type, 0.01 J7cm, monocrystalline Silicon wafers as pre- 
viously described in Refs. [26|, [27j]- The resulting porous 
substrates, covering an area of about 2 cm 2 , have average 
pore diameter d = 10—20 nm and average interpore spac- 
ing (defined as the distance between the centers of two 
consecutive pores) a — 20 — 50 nm. Nb ultrathin films 
were deposited on top of PS substrates in an UHV dc 
diode magnetron sputtering system with a base pressure 
in the low 10 -8 mbar regime following the same fabri- 
cation procedure described elsewhere prl 27]. Since the 
effect of the periodic template is reduced when the film 
thickness, d^b, exceeds the pore diameter [26(, the thick- 
est Nb films deposited for transport measurements have 
dNb = 12 nm. With this approach an array of intercon- 
nected Nb wires is guided by the PS template, the single 
wire width being the substrate pitch (w — 10 — 35 nm). 
Finally, the samples were patterned by standard optical 
lithography and lift-off procedure into bridges of width 
Wb = 10, 20 fim and length Lb = 100 /zm, to obtain a 
classical pseudo-4-point geometry (meaning two contacts, 
each used for a current and a voltage lead). The effect 
of the patterning is to reduce the number N of intercon- 
nected wires under study. 
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FIG. 1: Top-view scanning electron microscopy micrographs of (a) a free nanoporous Si substrate (a = 50 nm, d = 15 nm) and 
Nb nanowire networks deposited on PS with nominal thickness of (b) 3.5 nm, (c) 7 nm, (d) 15 nm. In the upper (lower) panels 
the images acquired at lower (higher) magnifications are reported. The white scale bar is 100 nm. 
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TABLE I: Summary of the calculated morphological param- 
eters at different film thickness (djvb). C indicates the grain 
size, while d is the pore diameter. Grains and pores average 
diameters are reported with an error of ±2 nm. 



Fig. Q] shows high resolution images of the surface 
of superconductive nanowire networks obtained from a 
template with a — 50 nm and d = 15 nm, performed 
by Field Emission Scanning Electron Microscopy (FE- 
SEM) (Sigma Gemini by Zeiss). The images of four dif- 
ferent samples are reported in the panels: (a) SO, a free 
nanoporous Si substate, (b, c, d) three Nb films (S3. 5, S7, 
and SI 5) with nominal thicknesses d^b = 3.5,7, and 15 
nm, respectively. The images were acquired at two dif- 
ferent magnifications (600 Kx and 1312 Kx) and clearly 
show the evolution of the morphology with the progres- 
sive Nb deposition. It emerges that the single nanowircs 
are polycrystalline with well shaped grains with increas- 
ing dimension of 10 - 20 nm for increasing dnb values. 
Images have been analyzed by SPIP [28], a software for 
image analysis, in order to evaluate the average diameter 
of pores, d, and grain size, These value are summarized 
in Table |U 

From the FESEM analyses, it follows that the op- 
timal Nb thickness to investigate the superconducting 
properties of the nanowire array should lie in the range 
d^b ~ 9 — 12 nm. This, in fact, assures that the wires 
are continuous, while preserving the presence of a well 
defined network. Moreover, at these reduced thicknesses 
the system could approach the ID limit. The character- 
istic parameters of these samples are listed in Table [TTJ 
Superconducting transition temperatures, T c , and criti- 
cal currents, I c , were resistively measured in a 4 He cryo- 



stat using a standard dc four-probe technique. Trans- 
port measurements provide also a probe of the quality 
and the homogeneity of the nanowires. It is worth un- 
derlining that on the same kind of unstructured samples 
an extensive electrical characterization was performed, 
which revealed that ultrathin Nb films deposited on PS 
substrates exhibit well established superconducting prop- 
erties [26l . [27J comparable with the ones of conventional 
plane films. It has been also demonstrated that the size 
and spacing of PS pores is well matched to the vortex 
lattice of Nb thin films, giving rise to matc hing effects at 
high magnetic fields and low temperature [26l [27|. The 
standard lift-off processing on a micrometer scale is not 
expected to induce drastic sample imperfections which 
could be responsible of a dramatic resistance versus tem- 
perature, R(T), broadening. Resistivity of the single 
wires is difficult to estimate accurately due to the par- 
ticular nanowires arrangements. However, on first con- 
servative approximation, it is reasonable to assume for 
Pn sl value comparable with the one obtained for plane 
Nb ultrathin films [29[ , which for sample deposited in the 
same sputtering system under the same conditions varies 
between pj^ = 50 and 35 pttcm for thickness in the ana- 
lyzed c/jv6 range [30(. Knowing the material constant pjyl 
for Niobium 3l|, namely pnI = 3.72 xl0~ 6 pttcm 2 , it 
results that the value of the low temperature mean free 
path is I w 1 nm in agreement with the values reported 
for ultrathin Nb nanowires [HJ . Moreover, being / -C £, 
it follows that the samples are in the dirty limit, an im- 
portant requirement for the occurrence of QPS processes, 
since the probability of a QPS event can be expressed as 

1/2 

P oc exp(—aT c ' ct/pn) (here a ss djyb-w is the wire cross 
sectional area) indicating that low temperature supercon- 
ductors with high resistivity are the best candidates for 
obser ving QPS. Compared to MoGe-based nanowire net- 
work [22| the detection of QPS is favored for Nb wires, 
since the required condition [l| fc^T < A(T), where A(T) 
is the temperature dependent superconducting gap, is re- 
alized at higher temperatures. Fig. [5] shows the resistive 
transitions for all the analyzed samples. Measurements 
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20 


10 
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9.5 


500 


10 3.39 


S10 


40 


30 


10 


17.3 


250 


10 3.52 


S12 


40 


30 


12 


19.0 


500 


20 3.57 



TABLE II: Samples characteristic dimensions: a indicates the 
interpore spacing, w the wire width, djvt the Nb thickness, <Ja 
the wire effective diameter, N the number of interconnected 
wires enclosed in the patterned bridge of width Wb, and T c 
the critical temperature of the samples. The numbers listed 
above are likely to be corrected considering a deviation from 
the corresponding mean interpore distance of the order of 
10% [2f|. The pore diameter, d, is 10 nm for all the measured 
samples. 



were performed using a constant bias current = 50 /iA. 
The first step in the curves is due to the transition of the 
electrodes [H, HU, therefore the values of T c , reported in 
Table HH were evaluated at the midpoint of the transi- 
tion occurring below the first one. Correspondingly, the 
normal state resistance, Rn, is defined as the sample re- 
sistance below the electrodes transition [H,[l6j]. The main 
feature of the R(T) curves is the nonzero resistance over 
a wide temperature range, which is strongly reminding 
of a ID behavior. It also deserves noticing that the R(T) 
curves do not show any ste ps o r humps, which can be a 
signature of inhomogeneity [ll|. Moreover, the moderate 
values of the normal state resistance (Rn <C Rq, where 
Rq is the quantum resistance Rq = h/Ae 2 s» 6.45 kfi) 
exclude that weak links constrictions or the film granular- 
ity can be responsible of the conduction [ll[ . In order to 
further demonstrate the good quality of the samples and 
to exclude the presence of tunneling barriers at the grain 
boundaries, it is useful to analyze both the shape of the 
voltage-current (V(I)) characteristics and the values of 
their critical current density, J c [l6| . Indeed at low tem- 
perature the shape of the V(I) characteristics (not shown 
here) are very sharp and present no steps or hysteresis. 
The values of J c are comparable with the ones measured 
on Nb perforated thin films 33 . The critical current of 
the network is, indeed, I c = (N + l)i c , where i c is the 
critical current of the single wire, which is related to its 
depairing current, id p , via i c — 2id p /3\/3. 34} It follows 
that from the measured I c values, the depairing current 
density of the single nanowirc at T — can be estimated 
to be Jdp — 1.5 x 10 10 A/m 2 , about one order of magni- 
tude smaller than the depairing current density evaluated 
in the framework of the model by Kupriyanov and Lu- 
kichev [H (KL) jg L = 3.5 x 10 11 A/m 2 , a discrepancy 
comparable to the one reported in the case of Nb perfo- 
rated thin films [33[ . Moreover, in the present case, this 
difference may be attributed, for instance, to the approx- 
imation performed calculating the values of N, namely to 
the complex morphology of the samples. Having clarified 
the issue of the sample homogeneity, it is possible to ana- 
lyze the R(T) curves in the framework of the theoretical 



models proposed for ID superconductors to describe both 
thermal and quantum phase slip processes. Indeed the 
systems under study consist of interconnected networks 
of ultrathin Nb nanowires with width (w = 10, 30 nm) 
comparable to the Nb superconducting coherence length. 
The latter was estimated from the temperature depen- 
dence of the perpendicular upper critical field, H C 2 ± (T) , 
obtained performing resistance versus field, R(H), mea- 
surements at fixed values of the temperature. Indeed 
it follows that the Ginzburg-Landau coherence length at 
zero temperature is £(0) = 10 nm. In presence of thermal 
activated phase slip processes the resistance is described 
by the equation [24j: 
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while QPS contribution may be expressed with exponen- 
tial accuracy as [H, HE HH : 
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where T* w T c is a crossover temperature between 
TAPS and QPS regim e, L is the nanowire length, and 
£(T) = Zo/y/l - T/T c and U(T) = U(0)(T c - T) 3 / 2 are 
the temperature dependent coherence length and phase 
slip activation energy, respectively. Since L is not well 
defined in these samples [22|], it is treated as fitting pa- 
rameter together with T c and, in the case of Eq. [T] ([2} , 
the parameter(s) U(0) (A and B), in the framework of 
the same approach used in Refs. |22|,|25j|. Thus, the total 
resistance of the samples can be expressed as 



R{T) = [R^ + [Rtaps + Rqpsr 1 ]- 



(3) 



The results of the theoretical analysis are reported in 
Fig. [2j For the sake of the clarity the best fitting curves 
are indicated by solid red lines. The data corresponding 
to sample S10 [Fig. HJa)] present a negative curvature 
consistent with thermal activation process, namely the 
lnR(T) curve does not level off at lower T as character- 
istic of quantum tunneling and they were satisfactorily 
fitted retaining only the TAPS term (Eq. [J) in Eq. E) 
Indeed, sample S10 is characterized by a slightly lower 
i?7V value and a slightly higher T c value compared to S9 
and S12. In fact, from Fig. [U[a) it can be inferred that 
including also the QPS term (Eq. ^ in Eq. [3] produces 
a much broader transition (dashed black line). 

The results of the theoretical interpretation according 
to Eq. |3]for samples S9 and S12 are reported in Fig. [^b) 
by thick red lines. As far as the sample S12 is concerned, 
the best fit to the R(T) curve obtained with Eq. [3] nicely 
follows the data, which, on the other hand strongly devi- 
ate from the solely TAPS dependence obtained from the 
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T(K) 

FIG. 2: (Color online). Resistive transitions, R(T), of dif- 
ferent Nb nanowire networks: (a) squares represent the ex- 
perimental data referring to sample S10, the solid red line is 
the theoretical curve obtained including only the TAPS con- 
tribution in Eq. [3 while dashed black line is the result of the 
fitting procedure including both TAPS and QPS contributions; 
(b) experimental data of samples S9 (circles) and S12 (trian- 
gles) are shown together with the curves obtained from Eq. [3| 
including only the TAPS term (dashed black lines) and both 
TAPS and QPS terms (solid red lines). 



same equation disregarding the QPS term (dashed black 
line). These results suggest that the pronounced resis- 
tance tail exhibited by sample S12 could indeed originate 
from QPS processes. On the other hand the R(T) de- 
pendence of bridge S9 arising from both TAPS and QPS 
(Eq. |3j is in agreement with the experimental points only 
down to T i=a 2.75 K. Below this temperature the curva- 
ture of the lnR(T) curve changes again with a slope more 
consistent with thermal activation. However, the discrep- 
ancy between the TAPS dependence obtained from Eq. [3] 
considering only the thermal contribution (dashed black 
line) and the experimental data is extremely pronounced. 
It is worth commenting the values of the parameters ex- 
tracted from the fitting procedure. In particular the val- 
ues of the critical temperature, namely T c = 3.5,4.0,3.5 
K, estimated for samples S9, S10, and S12, respectively, 
are close to the measured ones (see Table [TTJ) . The best 
fitting procedures performed with the full Eq. [3] was ob- 
tained using B sa 2 — 40, in agreement with the values 
reported in the literature (25|, but at the price of con- 



sidering quite high values of A sa 1000. This discrepancy 
can be possibly attributed to the exponential approxima- 
tion assumed in treating the QPS contribution. Finally 
for sample S10 a value of U(0) = 2.5 meV was extracted 
[H. From the relation U(0) sa 0.83LR Q T C /£(0)R N valid 
for a single nanowire in the framework of the approxi- 
mate phenomenological formula of Arrhenius-Little [25| . 
following Ref. [22[, it is possible to estimate L sa 312 ran. 

This analysis reveals that, despite the models were 
derived for an individual wire, they can reproduce the 
transition of networks with not well defined length and, 
moreover, with finite widths and activation energies dis- 
tribution. Clearly these spreads are not taken into ac- 
count in these models, as well as the distribution of the 
strengths of the contacts between the nanowires in the 
network. Indeed poor contacts, which may be present 
in these systems for instance at grain boundaries, can 
behave such weak links, the presence of which increases 
the probability of TAPS and QPS [H|. These param- 
eters fluctuations, which are due to the peculiar growth 
technique of the network, are expected to affect more the 
samples with lower values of cr 1 / 2 . In this sense the R(T) 
behavior of sample S12, which is expected to be more ho- 
mogeneous compared to sample S9, is better reproduced 
by Eq. [3] In this respect, however, it is worth comment- 
ing that, since in arrays of small Josephson contacts the 
thermal activation of the Josephson weak links may also 
affect the shape of the R(T) curves, the experimental 
data were also analyzed in the framework of the theory 
of Ambegaokar and Halperin 36] within the formulation 
given in Ref. @. The procedure revealed an extremely 
poor agreement between the R(T) dependencies and the 
model for any reasonable physical value of the fitting pa- 
rameters, confirming that weak links do not dominate 
the transport in these systems. Finally, to shed a light on 
the physics revealed by these systems two complementary 
experiments could be performed. On one hand electron 
beam lithography could be used to fabricate much nar- 
rower bridges, obtaining better controlled arrays with a 
drastically reduced number of interconnected wires and, 
consequently, a narrower distribution of widths and con- 
tact strength. Moreover the analysis of the microwave 
response of the arrays could help to determine the dom- 
inant dissipation process, namely TAPS or QPS [25j]. 



In summary, both thermal and quantum fluctuations 
were revealed from R(T) measurements in superconduct- 
ing Nb nanowire networks patterned on PS substrates. 
The innovation of the analyzed systems lies both in the 
widely accessible fabrication technique, which rely on 
rigid templates, and, most importantly, in the multiple- 
connectivity of the wires. 
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